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Abstract 

o\ 

f^ ' Throughout more than two millennia many formulas have been obtained, some of them 

f^ I beautiful, to calculate the number tt. Among them, we can find series, infinite products, 

04 ' expansions as continued fractions and expansions using radicals. Some expressions which 

^H I are (amazingly) related to tt have been evaluated. In addition, a continual battle has been 

^H' waged just to break the records computing digits of this number; records have been set 

■^H I using rapidly converging series, ultra fast algorithms and really surprising ones, calculating 

^sq ■ isolated digits. The development of powerful computers has played a fundamental role in 
these achievements of calculus. 

o' . 

}nr] . 1 First formula: Archimedes' algorithm 






For a period of approximately 1800 years, Archimedes' algorithm, proved around 250 BC, was 
the most efficient way to calculate vr. The idea consisted of considering a circle with unit diameter 
with regular circumscribed and inscribed polygons of 3 • 2" sides. Let a„ and 6„ denote their 
perimeters, respectively. He proved the following relations using only geometrical reasoning: 

m 

^ ' ^ r 2a„6„ / ; 

(N ! ai = 2V3, Oi = 3, a„+i = — — , bn+i = \^ an+iK- (1) 

t-^ . Qn + bn 

00 

Q . Obviously, 6„ < tt < a„ and both sequences a^ and hn converge to vr. It is an algorithm which 

I I \ nowadays can be easily proved using elementary trigonometry. Indeed, if /cn = 3-2", we can 

O ' write 

^■^ I -TT -TT 

OO . an = K tan — , 5„ = kn sin — . 

\^ • Kn Kn 

Then 

k> 2an6„ _ tan ^ sin ^ _ tt _ 

rS ' n~ ~ 2«nT F"~ '■ ^ — 2Kn tan — — — a^+i. 

\Jt • On + On tan T- + sm T^ 2kn 

On the other hand. 



/ 7~ /r.7 "^ r ■ "^ 1^7 vr.Tr vr , 

y/an+ibn = J^kn tan — - • kn sm — = 2A;„ Wtan — - sm — - cos — - = b„ 



ZKn Kn V '^Kn ^Kn Zk^-^ 



-1- 



Convergence in this algorithm is linear; after 5 steps we get 3 digits of vr. To see this, we use 

2 

the inequalities tan x < x and 1 — cos x < ^ to get 



2 
O-n — bn = kn tan -;— ( 1 — COS -;— ) < :rT7r = T^ — TTT- < 



vr / vr \ vr'^ vr'^ vr^ 



kn\ knJ 2kl 18 • 22" 18 • 1000"/^ ' 

1 



With n = 7, this algorithm gives the estimate 3.1415 < vr < 3.1417. 

It is enhghtening to point out that during the age of Archimedes neither decimal notation 
nor any other positional notation was known, so he stated the result he obtained in terms of 
polygons with 96 sides (n = 5) using fractions: 

10 1 

3+^<vr<3+-. 

This approximation determines two-place accuracy; this means vr ~ 3.14. 

2 No great news until 1800 years later 

For a period of approximately 1800 years all formulas that were obtained for tt were based on 
the idea of approximating a circle using regular polygons, and they took the shape of iterative 
algorithms. In 1593 Viete obtained a different expression by means of a pretty expansion with 
infinitely many nested square roots: 




11/11/1 
vr V2 V2 ' 2V2 ^ 2 "^ 2 V 2 "^ 2 V 2 " " ' ^ -* 

As with all the previous methods, Viete's is based on geometry and consists of considering a 
circle whose diameter is 2 and approximating its area, that is, vr, by the area of an inscribed 
square, the area of an inscribed octagon, the area of an inscribed 16-sided polygon, and so on. 
In this way Viete found vr accurately to nine places: 3.141592653. 

Around 1650, analytical methods of Descartes and Fermat, who were forerunners of the 
development of Calculus which would be carried out soon, started replacing geometrical methods. 
One of the formulas for vr at this stage is the famous infinite product due to Wallis (1655): 

2_13355779 

^ ~ 2 " 2 " 4 " 4 " 6 " 6 " 8 " S"""' ^ -* 

In 1999, Osier proved a formula which includes as particular cases both Viete's formula and 
Wallis' infinite product. We will see this in section 10. 

3 The early stages of the Calculus: New Formulas for tt 

Between 1665 and 1680, Newton and Leibniz independently developed the first basic notions of 
Calculus, which supplied us with a powerful tool not only for mathematics but also for physics 
and other sciences. These methods, as it could not be otherwise, also had repercussions in 
obtaining new formulas for vr. For example, one obtained by Newton is 



TT 



^ + 2-^y ffl (4) 

4 A^^{n + l){2n + 5)16"'' ^ ' 



by which he approximated vr to 15 places in 1665. Newton reached this formula by noticing that 
the area which is bounded by a circle with the equation 



1\^ 2_1 



2:- 2 ) +y 

and the straight lines x = j and y = is equal to the difference between the area of a circular 
sector that takes up an angle of 60° and its corresponding triangle. Therefore 

rl/4 /■1/4 /■1/4 ^ ^ 

/ V dx = V X — x'^dx = X ' (1 — x) ' dx = ^ . 

Jo ^ Jo Jo ^ ' 24 32 

Newton's formula is obtained from here by expanding (1 — x)^'^ using the binomial formula, 
which he had proved, and then integrating term by term. Formulas such as (jlj are nowadays 
known as binomial sums. 

A series similar to Newton's one but simpler is 

1 = \- ^n' (5) 

3 ^ (2n + l)16"' ^ ^ 

n=l ^ ' 



which we get by substituting x = \ into the series expansion (see |7j, p. 386) 

arcsin 2x 
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00 ('2n\ 



2x 



E \nl 
-, 2n+l 



x^-. 



n=l 

Another interesting example of a binomial sum is the following obtained by Comtet in 1974: 



i 



3240 ^ 1 



n=l \n I 



which can be deduced from the identity [H] 

/■7r/3 



1 l"^i'^ / ^\ 



In 1671 Gregory integrating the following geometric series term by term: 

l-x^ + x'^-x^^ = - ^, \x\<l, 

1 + x^ 

got the formula (which has his name) 

Q c; 7 

XXX II. 

arctan x = x 1 1----, x<l, 

3 5 7 ' II-' 

which would become known as the power series expansion of the function arctan x. Substituting 
X = 1, Leibniz got in 1674 the series for vr 



It is strange that Gregory did not do it himself, probably because he was not interested in slowly 
convergent series. Sharp put x = v3/3 in Gregory's formula and obtained the series 



^ _ \/3 / _ ^ _^ l__ \ 

6~3V 3-3'^5-32 7-33 "^"V 



(8) 



which allowed him to achieve the record of 71 digits for vr in 1699. It is surprising that Gregory's 
formula and Sharp's particular case have been found in the Sanskrit texts that the Indian learned 
person Nilakantha Somayaji had written 200 years before. 

In 1665 Brouncker got the pretty expansion as a continued fraction 

4 12 

- = 1 + ^2 • (9) 

2 + ^ 



2 + 



52 



2 + 



2+'-. 
About 100 years later, Euler proved it in a clever way. He starts by proving the equality 

ai + aia2 + 010203 + 01O2O3O4 • • • = ^- , 

1 



1 + 02 



1 , '^4 
1 + 03 



1 + 04 - ■ ■ . 

which motivates him to rewrite Gregory's series as 

(-x^\ (-x^\ /-3x2\ /-x2\ /-3x2\ /-5x2\ 

and then he gets Brouncker's formula by substituting x = 1 and using the continued fraction. 
A recent example of expanding vr as a continued fraction is 

^ = 3 + ^2 ' (10) 

6 + 



6 + 



52 



6 + 



6 + 



92 



6 + 



obtained by Lange [22] in 1999. 



4 Machin type series 

In 1706, John Machin proved the following formula (which is known today as Machin's formula): 

vr , 1 1 ,_, 

— = 4 arctan arctan . ( 11 ) 

4 5 239 ^ ^ 

Using Gregory's formula he rewrote it as a series and in this way he got the first 100 digits of 
vr. The following identities are similar, but they lead to more slowly convergent series: 

TT 1 1 vr 1 1 

— = arctan — h arctan — , — = 2 arctan — h arctan — . (12) 

4 2 3 4 3 7 ^ ^ 

The proof of a known Machin type series expression consists of an ordinary checking of an 
identity between elementary operations with complex numbers. We use 

^ ' =2 + 2i 
239 + i 

to prove Machin's formula; it is enough to make the arguments equal. 

Gauss (1777-1855) also took part in the search for rapidly convergent series for vr. He proved 
the following equality: 

vr 11 1 

— = 12 arctan h 8 arctan 5 arctan . (13) 

4 18 57 239 ^ ' 

It is like the preceding ones, but it involves one more term. 

In 2002 Kanada and his team, using a certain type of trick to optimize computer calculations, 
got a record computing vr by means of the series 

TT 1 1 1 1 

— = 12 arctan h 32 arctan 5 arctan h 12 arctan , ( 14) 

4 49 57 239 110443' ^ ' 

— = 44 arctan h 7 arctan 12 arctan \- 24 arctan . (15) 

4 57 239 682 12943 ^ ' 

The first one was obtained by Takano in 1982 (he is a high school teacher and a song composer). 

The second one was obtained by Stormer in 1896. Kanada's record consisted of more than a 

billion of digits of vr. 

5 Formulas by the Swiss raathematician L. Euler 

Great mathematicians, Leibniz and the Bernoulli brothers Jacob (1654-1705) and Johann (1667- 
1748) among them, had attempted to evaluate (with no success) the infinite sum 

11111 
Finally, Euler proved in 1736 that its value is 7r^/6. To do it he started with the series 

2 4 fi 

smx X X x" 



^ 3! ^ 5! 7! ^ 



and the equation 



2 4 fi 

X X X 



^-3r+5!-7!+--- = °- 



The solutions of this equation should be the solutions of sinx = 0, except x = 0, that is: ivr, 
±27r, zLStt, .... He imagined the infinite series as a polynomial and he broke it down into factors 

x^ x^ x^ ( x^\ ( x^\( x^\ 

He finally got the proof by setting the x^ coefficients equal. He also proved many similar 
formulas. Among them we can find 

^(2n+f)2 8' ^ n4 90' ^, (2n + 1)3 32 ^ ' 

n=l ^ ' n=l n=0 ^ ' 

and the not-so-similar 

. 1 1 1 

oo in — I — I 4 

V 2 3 n _ ^ ns^ 

^ ^3 ^- (18) 

n=l 



Nowadays formulas like the last one are known as Euler sums [15]. The following one has an 
interesting history 

/ 1 1 1\2 

V A 2_^ VJ_ - iZ^ (19) 

^ V? 360 ' ^ ^ 

n=l 

because it was observed numerically by the student Au Yeung. J. Borwein, his professor, was 
skeptical of the result but after checking it with greater precision he got to prove it. Now, 
we know that the formula had been proved long before by H. Sandham [30]. However, the 
rediscovery of this forgotten formula has had an important effect because it has motivated the 
organized study of these kind of formulas. 

While studying the divergence of the series of reciprocals of the prime numbers, Euler con- 
sidered geometric progressions: 

1111 1 



p2 pi p& pS -^ ^ ' 

and he observed that by multiplying over all primes p we get 

n=l P ^ 

which links the number tt with prime numbers. 

From 1748 he initiated the study of the lemniscate and discovered the formula 

dt /■! t^dt TT 

( Ll 1 



/o \/r^ h \/r^ 4' 

which establish a remarkable relationship between two elliptic integrals. 



6 An outstanding forgotten algorithm: C. F. Gauss' Formula 

When he was only 14 years old, Gauss considered the following iterative process: 



oo = a, bo = b, a„+i = , bn+i = yo-ubn 

and he defined the arithmetic-geometric mean of a and b as the common limit 

M(a,b) = lim a„, = lim 6„. 

n— >oo n— >oo 

In 1800, his investigations concerning the arithmetic-geometric mean lead him to the wonderful 
formula (see [l], p. 94-102) 

E2^(4-fti) = l-^ (22) 

— TT 

where ao = 1, 6o = l/v2 and M is the arithmetic-geometric mean of ag and bo. 
We sketch a proof of it using the elliptic integrals 

^/2 de ^ r/^ sin^ede 



J-n — / I ; J-in 



a^ cos^ 6 + bl^ sin^ 6 "'o J a^ cos^ 6 + b'^ sin^ i 



The first step consist in proving that /„ = /«+! by means of clever substitutions of the integration 
variable. This result has an important consequence, namely: 

TT 

Jo = In = hm /„ = — — . 

ra^oo 2M 

The second step consist in deriving, again by clever substitutions, the relation 

clLn - 2c^+iL„+i = -clln, where cl=al- 6^. 
Multiplying by 2" and summing for n > 0, we obtain 

-, oo 

^ n=0 

The last step is the evaluation /o(/o ~ -^o) = '^/2, which follows easily from Euler's formula (I2ip . 
For no clear reasons. Gauss did not make use of his formula to calculate tt. In 1970, it was 
rediscovered independently by Eugene Salamin and Richard Brent |26j, [TU] who expressed it 
appropriately for iterative calculation: 

hm = TT. 

n— >oo _ ". 

l-Y.'i\al-bl) 

k=0 

It can be proved that the convergence to vr is quadratic, that is, from n to n -|- 1 the number 
of correct digits is doubled. Brent and Salamin programmed this algorithm on a computer. It 
let them reach in 1976 a record approximating vr to 3 million digits, while the record until that 
moment had been 1 million digits, achieved using Machin-type formulas. 



7 Series for I/tt from an incredible Indian mathematical genius: 
S. Ramanujan 

The perspicacity and ability of Ramanujan allowed him to discover, among many other things, 
17 extraordinary series for l/vr completely different from those known until that moment. These 
new series were published in 1914 in his celebrated paper |25j . In it he considers the complete 
elliptic integral of the first kind 



Jl-B^-^-- V2y V2.4 



' sm 

The number k is called the modulus and k' = \/l — k"^ is called the complementary modulus. 
He then gives his definition of what ordinary modular equations are: They express the algebraic 
relations which hold between k and I when, for a given rational r, we have rV /L = K' /K, or 
Q^ = q, where q = e~'" ' , Q = e~^ ' and X, K' ^ L and L' denote complete elliptic integrals 
with the moduli k, k' , I and I' . Supposing k = I' and k' = I, we have q = e~'^^^. Developing 
this theory [23], [5] he proved the formulas 

n=0 n=0 

and 

n=0 \ / 

He says that the other 14 examples he gives belong to three alternative theories: the quartic, 
the cubic and the sextic. He defines the corresponding functions K, but it is a pity that he does 
not develop these theories. One of his examples in the quartic theory is 

^ (4n)! 1 , , 9801^/2 , , 

Y^ ^ 26390n + 1103) = — -^— . (25) 

n=0 

In the cubic theory he gives, for example, the series 

^ (2n)!(3n)! 1 , , 27 

y ^^ — ^T^^^ 15n + 2) = — . (26) 

n=0 

And in the sextic 

^ (6n)! 1 , , 5\/l5 , , 

y , ^ ,/., (lln + l) = ^^ — . (27) 

n=0 ^ ' 



Series ()25p is incredibly fast: each term contributes 8 digits of vr. In 1985, before a proof of this 
formula was known, W. Gosper programmed it on a computer and got a record of 17 million 
digits of vr. 



The Borwein brothers, Jonathan and Peter, were the first to give complete and rigorous 
proofs, pubhshed in 1987, of the 17 Ramanujan series [7]. At about the same time the Chud- 
novsky brothers, David and Gregory, developed in |13j a theory to prove the following amazing 
series formula, which was not discovered by Ramanujan: 

^ (6n)l (-1)'^ , , \/6403203 , ^ 

Y t^t/tt^^ ^ 545140134n + 13591409) = . (28) 

^ 3n!n!^ 6403203'^ ^ ^ 127r ^ ^ 

n=0 ^ ' 

Each term contributes 14 digits of vr and they used it to beat records of the calculation of vr 
in 1989, 1991 and 1994, the last one with 4044 million digits. In addition, they proved that 
there formula (j28p . which corresponds to r = 163, has the fastest convergence rate among all 
those Ramanujan series that involve a rational raised to the power n. They arrived at this 
conclusion by noticing that d = 163 is the largest positive integer such that the field Q{\/—d) 
has the unique factorization property [IS]- In the year 2001, H. H. Chan, W. C. Liaw and V. 
Tan found identities [llj which allowed them to prove the following formula, which also had not 
been discovered by Ramanujan: 



(3n)!(2n)!(-l)",^^^^^^ , oov^ _ 1500^3 

ra=0 



s^^^^yfe*"'^'"^^^^'-^^- p'> 



Very recently N. Baruah and B. Berndt in [3j and [4j following more closely the initial ideas of 
Ramanujan as presented in Section 13 of the celebrated paper [25j got a much simpler treatment 
of the theory, and proved many new series of this type. An excellent survey on Ramanujan's 
series is [5J. 

8 The Borwein brothers' quartic algorithm 

The fourth root of the elliptic lambda function and the elliptic alpha function [7], that is, 



,{qy ^^^ vr eliq) 



^-(^''-)|i 



where 



n=oo n=oo 



n=— 00 n=— 00 n=— 00 

are the Jacobi elliptic theta functions, satisfy the following wonderful properties: 



^(^4. _ 1-^1^^ 



l+^l-s4(g)' 
a{q^) = aiq)[l + siq^T + ^^siq^)[l + siq^) + s\q% 



IT 

These identities connect the function values at q with the values at g^ , so they allow us to obtain 
an amazing Borwein family of quartic algorithms ^. For example, if we let 



q = e ""' , Sn = s{q), tn = a{q), 



then the values of sq and to are algebraic, and we get the following quartic algorithm: 



Sn+l 



i-n 



l + 0r^l: 



¥ 



(30) 



f^n+l 



tn{l + Sn+l)^ - 22"+3s„+i(l + Sn+1 + S^ 



where tn tends to l/vr, and the number of accurate digits is quadrupled in each iteration. In the 
same way, if we let 

q = e~'^"^", Sn = s{q), tn = a{q), 

we obtain a quartic algorithm with initial values sq = 1/ v2 and to = 1/2- A simple proof of this 
one is given in [21j using only Gauss' formula (j22p and elementary algebra. However, although 
this proof is much simpler than the Borweins', there method is more general. In addition, J. 
Borwein and P. Borwein have obtained other quadratic, cubic, etc algorithms. However, it seems 
that the quartic algorithms are the most efficient. 

We might think that vr calculation records would always be broken using quadratic, quartic, 
etc. algorithms. However, that is not true, because such algorithms do not auto-correct their 
errors and they require doing all the operations to the final desired accuracy. Other ways have 
been found to improve the accuracy of calculation of series without affecting the final result, 
allowing them to beat records even though they have the drawback of being only linear. 

9 The jumbled calculation of tt digits 

In 1995 Peter Borwein and Simon Plouffe realized that the series 

_ 1 1 1 1 _ '^ 1 

n=l 

could be used to compute one of the binary digits of In 2 without calculating the previous digits, 
by means of the following ingenious formula, where {x} means the fractional part of x: 

( ~i I ( '^ 2''~"1 °° 1 1 f f ^ 2''~"modnl °° 1 

|2'=in2} = ||x:--|+^|:^^-^^| = |l|: - j+^e^^^^^^ 

It lets us efficiently obtain the binary digit of In 2 in position k + 1. The main remark to carry 
out all the involved calculation is that the numerator of the first sum, that is, 2^~"' mod n, can be 
computed very quickly. For example, the calculation of 2^^ does not require 65 multiplications 
but only 7 if we do it as ((((((2^)^)^)^)^)^) • 2 and a very efficient algorithm is used to obtain for 
example 2^^ mod 7. It consists of reducing modulo 7 after each multiplication instead of doing 
them in a row. 

In 1989 the Borwein brothers claimed that obtaining just one binary digit of vr could not be 
easier than computing all the previous digits. However, after this discovery for In 2, it was not 
so clear and the team consisting of D. Bailey, P. Borwein and S. Plouffe began to search, using 

10 



the PSLQ algorithm, for a series for tt that would allow them to calculate isolatedly its base 2 
digits (see [2], chapter 8). 

The PSLQ algorithm is a numerical algorithm which operates to a prearranged accuracy. 
It links an output of n integers ai, 02, . . . , a„, where aj 7^ for any i, with an input of n real 
numbers xi, 2:2, • • • , Xn, and 

aixi + 022:2 H h anXn = 

is verified. 

Success arrived taking 



00 



' „?o(8^ + j')16"' ^ l'2'---'7, 



and trying to find, using PSLQ, a linear relation with integer coefficients using as input vr, xi, 
X2, X3, X4, X5, xq, xj. The output obtained was this sequence of integers: —1, 4, 0, 0, —2, —1, 
— 1, 0. That is, they found a series P] which they called BBP type (BBP are the initials of the 
team members' names) 

_ ^ 1 / 4 2 1 1 \ 

^ ~ ;^Q 16" Ura + l ~ 8n + 4 ~ 8n + 5 ~ 8n + 6y ^ -* 

that lets us calculate vr's hexadecimal digits (hence also its binary ones) in isolation. Later, a 
proof of this formula was obtained from the easy result 



Jo l-xs'"'^ io „t^o 2V\^,16«(8n 

Namely, we have 



^1/ 4 2 1 1 \ _ /■ V^ 4^/2 - 8x^ - AV2x^ - 8x^ 

;^g 16" V8n + 1 ~ 8n + 4 ~ 8n + 5 ~ 8n + 6/ ~ io 1 - x^ ^' 

and this last integral can be easily evaluated by substituting y = x-v/2, as follows: 

/■! 4 - 2j/3 -y^-y'' ^ [^ I6y - 16 

'Vo 1^^? '^ = 70 iy^-2)iy^-2y + 2/y = ^- 



As soon as they announced the amazing formula ()3ip . other investigators wanted to find 
other beautiful BBP type formulas. The easiest one for vr was found and proved by Viktor 
Adamchik and Stan Wagon (see [29] pp. 66 and 67) 

^ 4" V4n + 1 4n + 2 4n + 3y ^ ^ 

n=U 

It can be proved similarly to the previous one, just taking into account the identity 



l + x4 7o ^0' ' 2*^/2 ;^o 4"(4n + /c) ■ 

11 



In 1997, Fabrice Bellard proved the formula 



- = 4E 



(_l)n 



n=o4"(2n + l) 



j_^ i-iy 



32 



+ 



+ 



64 ^ 1024" V4n + 1 4n + 2 4n + 3 



(33) 



and used it to brake a calculation record finding that the billionth binary digit of vr is 1. His 
formula is now used in almost every computation of isolated binary digits of vr and the actual 
record belongs to C. Percival. 

There are also known BBP formulas for other values related to vr, for example. 



vr 



Vs 



9 ^ 1 



E 



16 



1 



32^64"V6n + l 6n + 2 6n + 4 6n + 5 



and 



q °° 1 

S ^-^ P.An 



16 



24 



6 



+ 



8 ^ 64" V (6n + 1)2 (6n + 2)2 (6n + 3)2 (6n + 4)2 (6n + 5)2 



(34) 



(35) 



10 Pretty infinite expansions for tt using products and radicals 

In 1999, Osier proved the following formula |23j : 



-=n 

vr ^^ 



1 1 



7: + ^\L + 



1 /I 



-^\2 2^2 2V2 '" 2V2 



+9V9-n 



00 9P+1 



2P+in - 1 2P+in + 1 



n=l 



2P+in 



2P+in ' 



(36) 



where n inside the square-root means the number of terms in the sum. This formula includes 
Viete's formula © (p = 00); Wallis' product ^ (p = 0) and new formulas. For example, with 
p = 2 he gets 

2 ,_._.,..„.,_„_. ^^^^ 



vr 



1/1 1 /I 7 9 15 17 23 25 
2V2'^2V2"8'8'l6'l6'24'24 



In order to reach his formula. Osier used the double angle sine formula p times in a row: 



sin X = 2 cos — sin — = 2 cos — cos -?r sin —^ = 2"^ cos — cos —^ 
2 2 2 22 22 2 22 



X X 

■ cos — sm — 
2P 2P 



and substituted i = ^ in Euler's formula 



sint 



n 1 



n=l 



7r2n2 



n 

n=l 



n-K — t niT + t 



nvr 



nvr 



Combining the preceding formulas he got 



sm X XX 

= cos — cos —^■■- cos — 

X 2 22 2P 



X -TT / 2Pnvr — x 2Pnvr + x 



n 

n=l 



2Pnvr 



2Pnvr 
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Finally, he used the elementary trigonometry relations 



a; /I 1 X 1 1 /I 1 



cos — = \ 1 COSX, COS —^ = \ — I — \ — I — cosx, 

2V22 ' 22V22V22 



and substituted x = ^. 



In 2002, J. Sondow obtained another pretty expression with products and radicals [27J: 

00 1 

^=11 [l(-i)Ho) . 2(-i)^(") . . . (n + 1)(-1)"+H:)] - . (38) 

n=0 

He started by taking logarithms in Wallis' product ([3]): 

ln^ = y(-l)'=ln^. 

Taking into account that J2n>k 2^1 = 1) '^s continue like this: 

inl = y(-i)*^in^ y -S) = y ^y (-i)fcR in^. 

fc>0 n>k n=0 k=0 \ / 

Then he substituted n — 1 for n, rewrote the logarithm as ln(A; + 2) — \n{k + 1) and the sum on 
k as the difference of two sums. He replaced k with fc — 1 in the first of these two sums. Then 
he applied the binomial relation (^Z^) + ("^ ) = (^) and by exponentiating the expression he 
got the result. 

11 Very recent results: Series for I/tt, I/tt^ and I/tt^ 

T. Sato, in a 2002 conference in Japan, showed the series 

yn^P-^— (20n + 10 - 3^5) = g^ , (39) 

where Un are the Apery numbers, defined by 

Soon afterwards this type of series was found for other types of numbers [12], [2H]. H. H. 
Chan and Y. Yang independently developed a theory explaining this kind of formula, which also 
includes as a particular case Ramanujan type series, and they named them Ramanujan-Sato 
type formulas. 



13 



A radically different explanation for some Ramanujan type formulas, one that does not use 
the theory of modular functions, is based on the WZ method, due to Herb Wilf and Doron 
Zeilberger [2l]. It lets one automatically prove identities of the shape 

oo 

22 G{n, k) = Constant 

n=0 

in case the function G(n, k) is hypergeometric in both of its variables, that is, if the quotients 
G{n + l,k)/G{n,k) and G{n,k + 1)/G{n,k) are rational functions of n and k. The program 
EKHAD, written by D. Zeilberger, does the hard work of finding a companion function F{n, k) 
to G(n, k) such that F(0, k) = 0, and which makes [F, G) a Wilf and Zeilberger (WZ) pair, that 
is, a pair characterized by the property 

F{n + l,k)- F{n, k) = G{n, k + 1) - G{n, k). 

In addition, Zeilberger has proved that if we define H{n,k) = F{n + l,n + k) + G{n,n + k), 
then we will obtain an identity associated to the first one: 

oo oo 

^ H{n, fe) = ^ G{n, k) = Constant. 

n=0 n=0 

It is important to emphasize that the WZ method is a symbolic method and not a numerical 
one like the PSLQ algorithm (see section 9), so the proofs obtained are completely rigorous. 
The problem is that it is a difficult task to discover suitable functions belonging to a WZ pair. 

Some identities discovered by the author [18] , among them 



oo . /2m3/2A:x2 

y , ^/"^ ^V (6n + 4fc + l) = -, (40) 



y 10 ^i. A (20^^ + 8n + 1 + 2Akn + SA:^ + Ak) = -^, (41) 

^^=0 \ k ) 

are suitable for having an automatic proof with EKHAD. Just to see how the WZ method works, 
we explain in detail the proof of the identity (j40p based on this method. Given the function 

-, /2m3/2fc\2 

Gin, k) = , ,u % (6^ + 4A: + 1), 

V k I 



EKHAD gives the companion 



1 o"o" 



F(n,k) = « , '"' 'V 8^ 

V k ) 

and together they form a WZ pair, that is, they satisfy 

F{n + l,k)- F{n, k) = G{n, k + 1) - G{n, k). 

14 



We sum over all n > 0, so that a telescoping cancellation is produced on the first member, and 
we get 

n=0 n=0 

This series converges uniformly, so we can write 

CO oo f2k\^ J 

lim y G{n, k) = Y 1™ ^(n, k) = lim G(0, k) = lim Hx-(4/c + 1) = - 

A:— >oo ■^"1 — 1 k—>oo k^oo k^oo 2 TT 

71=0 n=0 

and we obtain the result 

~ 4 

?i=0 

And, as explained before, we have the double identity 

oo oo ^ 

Y,Gin,k) = J2Hin,k) = -, 

n=0 n=0 ^ 

where H{n, k) = F{n + 1, n + A;) + G{n, n + k). Substituting A; = 0, we obtain the Ramanujan 
series (|23l) 

f (5!|!4(6„ + l, = lfe=|!^(42„ + 5) = i (42) 

/L. ^16 28" ^ ^4 ^ n!6 212" v ; ^ v ; 

n=0 n=0 

For the function G{n,k) in the series (141 p . EKHAD finds its companion 
nn, k) = ^\,„^s/y^y • 8n(2n + 4fc + 1), 

and then we can say that the proof of ()4ip is hidden inside this pair [F, G) . We have the double 
identity 



Y,Gin,k) = Y,Hin,k) = —. 

— — TT 

n=0 n=0 

Substituting /c = 0, we obtain the new formulas [TB], [TH], [20] 



^ (2n)!5(-l)", 2 ^ 1 ^ (2n)!5(-ir, 2 s 8 , , 

E ^^(20n^ + 8" + 1) = Y^ E ^ V^(820n2 + 180" + 13) = ^. (43) 

n=0 n=0 

Using the PSLQ algorithm, the author managed to find more series of this type for I/tt^ 
[TT] , for example, 

^ (6n)! (-1)" , 2 ^ 128^5 , , 

n=0 

whose proof is still unknown. In [T9] and [20] we found the same series as the result of a 
conjecture. 
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B. Gourevitch found a series for l/vr^ [17] using the PSLQ algorithm 

E ^i(168n^ + 76n^ + ^^n + 1) = f • (45) 

n=0 

It has not been proved yet. 

W. Zudihn proved the following new identities [31] by means of quadratic transformations 
of the identities ([13|) : 

^ (4n)! , ^ ,1 10^/5 , , 

^ Wn ,o/ \, (18?^ -10^-3 ,« ^, = ^-, (46) 



00 



V- (4n)! ,0 ,1 5^41^/41 , , 

E ^'^;^y2(^(1046529n2 + 227104n + 16032)^^^^^ = — ^, (47) 

where the Wn are numbers defined by 

3 

Wr. 



tp) f;:::>-". 



With these last formulas we conclude our historical review. 
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